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C. D. U. 532.517.2 


Movimento lamelar dum liquido viscoso 
em regime não permanente 


ror FREDERICO MACHADO 


Engenheiro Civil (1. S. T.) 


As equações gerais do movimento dum líquido viscoso são (cf. Lamb, 1945, p. 577) 


Du 1 dp 
EE 
[o P dx 
Dv 1 dp 
1 —=Y-—=+toyvk 
(1) à D ? dy , 
Dio 1 dp 
——“ =Z——— DA vdd o 
(Ds da 


com D/Dt=o0/dttud/dxt-vo/dy+t-wd/dz e vi=09/[dX+4d/dy+d/dz 


Aqui o é a densidade (massa específica), p a pressão, e v a viscosidade cinemática; u, v, w 
são as componentes da velocidade no sistema de coordenadas cartesianas ortogonais Oxyz; ve» 
consideram-se constantes em todo o líquido. 

O movimento deve satisfazer também à equação de continuidade (para líquidos incom- 
pressíveis). 


(2) Apa 


dx dy dz 


Consideremos um tubo circular uniforme de raio a e eixo vertical coincidente com Oz, ligando 
um depósito de largura infinita à profundidade z == — h a uma superfície livre também inde- 
finida 7 =0. 

Teremos X==Y= 0; e, supondo o escoamento lamelar, u==v=0. A primeira e a 
segunda equações de (1) dão dp/dx = dp/dy==0, o que mostra que p só depende de t e z. 
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Por sua vez a equação (2) reduz-se a dw/dz==0 e portanto w é só função de t, xr e y, 
ou melhor, por razões de simetria, função de te de r=>Vt +y*. 

Além disso é Z==—g e vamos fazer p = — tgz + Ap, sendo 4p a diferença de pressão 
relativa à pressão hidrostática —-gz. A terceira equação de (1) vem 


dt o 


y, 
dw e dSp F vvtio 


EE dz 


p e Ôp são funções de te z, mas a equação anterior mostra que dAp/dz é só função de t 


EA . 


(porque w é independente de z). Será então 4p == Az +4-B (sendo 4 e B funções só de t) que 


deverá satisfazer às condições Ap ==0 para z=-0, e dp = Ap, para z=—h. Podemos escre- 
ver Ap=— (z/h)Ap, e — (1/g) ddo/ dz = Apulçh = E (t). 
Fazendo 


f 
w=w— [E (1 dt 
to 


(onde to é uma constante a determinar oportunamente) e passando a coordenadas cilíndricas, com 
w dependendo apenas de t e r, vem finalmente 


OW (- IV 1 0W 
3 —— = y - O ea 
(3) dr? + r dr 


com a condição nos limites 


f 
w== 0 eportanto IV = — fe (+) dt para r=a 
to 


(o que equivale a supor que o líquido não escorrega sobre as paredes do tubo). 
Uma solução de (3) é (cf. casos análogos indicados por Byerly, 1959, p. 12 e 226) 


(4) W=qe Jo (rVca) 


onde ci e c; são constantes arbitrárias, e Jo (rVc2) é a função de Bessel de primeira espécie e de 
ordem zero 


6 


— 2 Gr 
O ad E SR 


* + x 


Para satisfazer à condição nos limites, convém transformar a expressão (4). Suponhamos que 
F (t) é uma função de variação limitada, e que pode representar-se pela série de Fourier (cf. por 
exemplo Jordan, 1883, p. 228 e seg.) 


ps nat 


(5) El)== 3 À, sen- 
Ed 


n==1 


com 


T - 
(6) A= [E a u 
E Já T 
onde T é uma constante. 
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Será também 


Ê oo - 
(7) feoa=—— Rd 
=] MN T 


to Tn 


onde escolhemos to de modo a anular-se, na série, o termo independente de t. 
Façamos agora em (4) q==cetivitr e = +Hizx com a==n7/vT. Virá sucessivamente 


W' =c[cosvz (t—t,)—isen va(t—t)]), (rVia) 
W'=clcosva(t—t)-+isen v:(t—t)]], (rV-ia) 
Somando temos 
(8) W=eW+W= cn [cos va (t— ti) L(rVa) — sen va (t— t,) M (rV2)] 
VEL (avo) + IM (ava) 


sendo 


o pá at y8 
22.42 Maes "T 


L (rV2) = = [J, (VT +, (rVDT)]=1 — 


a3 48 


r 
2? “OA e 


ay 


M (rVa) = > (Via) — Jo (VD 0)] = 


c=0/2VI L(aVo)]? + [IM (aVa)) 


Estas expressões são reais, apesar das funções de Bessel terem argumentos imaginários. 
Se tr for determinado pela equação 


M (aV'«) 


t va t, == = 7 
. L (aVa) 


a expressão (8) dará, para r=a, 


nat 


! 


(Wc) == Cu COS VAL == é, COB 


o que permite escrever 
( Wa), = = 


com wm=(T/7) A,/jn e A, dado por (6). 
Então se fizermos 


Ar — Cos vz (E— &,) L(rVx) — sen va (tt) M (rVa) 
n 


(9) W = z W, == E = 
VIL (aV5)]? + IM (aV'a)]? 


oo 
2 


1 


com «=n7/vT e t=(1/v2) arctg [M (a Va) /L(aVa)], teremos, para r=a, 


A expressão (9) é portanto a solução da equação (3), satisfazendo à condição nos limites. 
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A velocidade do líquido será 


(10) o 7 | EA 


to 


e o caudal correspondente 


= 27 “wrár. 
(11) Q f 


Este problema de escoamento lamelar dum líquido viscoso, através dum tubo rectilíneo de 
secção circular uniforme, fica assim completamente resolvido. 


A solução pode porém simplificar-se se for suficientemente pequena a derivada ddpu/dt. 
Atendendo a (5), temos 


Â ss % 
Pa = F(t)= 5 As sen "E 
ch n=] T 
e, segundo (6). 
E 2 


nº 
So Aps sen ae dt 


2 E n=t 
p= (+) sen Ec 


expressão que, segundo a teoria das séries de Fourier, tende para zero quando n tende para infinito- 
Integrando por partes obtemos 


Apu 
An =— K Apn )i=o0 E ( Sp; )t=T cos nx he E ais cos e a 


ai 


Como por hipótese Spy tem variação limitada, se d Apu/dt for suficientemente pequena pode- 
mos então considerar valores grandes de T, sem que os A, aumentem apreciâvelmente, de tal 
forma que (nx/vT)r? ou (nx=/vT)a? sejam pequenos em face da unidade (podendo despre- 
zar-se as potências de ordem igual ou superior à segunda) para aqueles valores de n inferiores ao 
limite acima do qual os A,/n são suficientemente pequenos para se poderem desprezar os termos 
correspondentes da série (9) onde os outros são todos limitados. 

Teremos, nesta hipótese, 


n= a 


sen val;=—. AR cos vatln= 1, 
vT 
L(rVa)=1, tis. 
vT 4 
SET = = VR 
VEL (aVa)]? + [M (aVa)]? = 
e virá 
e FR As A. Du pt O r sã 
W=— 2 cos —— + CTT > Ansen” im [ Ergo É e 
Fu=] n T 4v u=1 4v “ch 


em consequência de (5) e (7). 
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CONSTRUÇÕES METALO-MECÂNICAS MAGUE, L.?* 


(S. E. MONIZ DA MAIA & VAZ GUEDES, L.º*) 
Alverca — Portugal 


Aspecto do invólucro espiral de uma das duas tarbinas «Francis» de 54.000 CV destinadas à Central 
de Pisões (HICA) construídas por Construções Metalomecânicas MAGUE S.A.R.L. segundo 
projecto e licença dos Ateliers de Constructions Mécaniques de Vevey S.A, 


PRONECTO E FABRICO DE 


PONTES ROLANTES 
GUINDASTES 
GRUAS. DERRICKS E GUINCHOS 
APARELHOS DE ELEVAÇÃO ESPECIAIS 
CONSTRUÇÕES METÁLICAS PESADAS 
TECNICA X 111 


- 


vnnmAs 


m 
3 
m 
=) 
e 


PARA 
ARMAZENAGEM 
AMONIACO 
AÇOS VAZADOS e FERROS FUNDIDOS 
MECÂNICA MÉDIA E PESADA | 
"CONSTRUÇÕES METÁLICAS 
comportas 
turbinas P 


condutas forçadas 
reservatórios sob pressão 
permutadores de calor 


REDUTORES DE VELOCIDADE 


3% 


DIVISÃO METALO-MECÂNICA 


COMPANHIA UNIAO FABRIL 


AV. 24 DE JULHO, 170 — LISBOA 2 — TELEFONES 670421 E 672061 — TELEGRAMAS: METALOFABRIL 
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Hu 


Substituindo em (10) e (11), obtém-se 


a* — 1º Ap; 
D == os! 
4v oh 
e 
x a! 
ai Apa 
8vgh 


que é a fórmula de Poiseuille (cf. Lamb, 1945, p. 585). 


Vemos que esta fórmula é verdadeira para Ap, constante (regime permanente), ou ainda 
aproximadamente, quando a pressão tiver variações suficientemente lentas (e puder desenvolver-se 
numa série da forma (5) ou análoga). O problema tem muito interesse para o estudo da ascensão 


da lava nas chaminés vulcânicas (cf. Machado, 1962). 


Agradecimento 


Ao Dr. F.R. Dias Agudo o autor agradece a leitura do manuscrito do presente trabalho. 


SUMMARY 


In this paper a solution is obtained for the variable flow of a viscous liquid through 
a vertical circular pipe (intended to represent a model of volcanic vent). When the pressure 
varies very slowly, the expression for the flux reduces to Poiseuille's law, as it should be 
expected. 
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APLICAÇÃO DO MÉTODO DAS ANALOGIAS ELÉCTRICAS AO ESTUDO 
DE ALGUNS PROBLEMAS RELATIVOS AO CAMPO MAGNETICO DAS 
MAQUINAS ELÉCTRICAS 

PELO ENG.º ELECT. (1.5.7.) CARLOS DE MEDEIROS PORTELA 


Assistente do 1.95. T. 
Continuação) 


5 Campo magnético na zona do entreferro, dentes e cavas 


5.1 Decomposição do campo magnético na zona do entreferro, dentes e cavas, para efeito da sua determina- 
ção experimental 


O campo magnético que mais afecta o comportamento electromagnético das máquinas eléctricas 
e cujo estudo é mais difícil, é o campo na região do entreferro, cavas e dentes, em que a influência 
da saturação do ferro, especialmente nos dentes, é decisiva. Este campo pode, à excepção de algumas 
pequenas zonas, ser tratado como campo a duas dimensões. 

Seja por exemplo o problema da determinação das reactâncias do induzido de um alternador, 
determinação que é possível reduzir à determinação do campo associado às correntes que percorrem 
o indutor e os diversos condutores do induzido. 

O campo numa zona do entreferro, cavas e dentes, associado a correntes no induzido não 
muito próximas, pode normalmente ser tratado como o campo associado às correntes indutoras, 
considerando a soma das forças magnetomotrizes do indutor e dessa parte do induzido. 

O campo associado às correntes que percorrem uma cava, na sua vizinhança imediata, não 
pode, porém, ser tratado da mesma forma, por ser totalmente diverso do campo associado às cor- 
rentes do indutor. 

O campo magnético na zona do entreferro, dentes e cavas, pode ser calculado da seguinte forma: 


a) Considerar separadamente o campo associado a cada dente e metade das cavas vizinhas, 
tendo em conta a corrente do indutor e do induzido em todas as cavas da máquina à excepção 
das cavas vizinhas do dente. Este campo pode ser determinado com bastante aproximação, 
considerando a saturação do ferro, a partir da relação entre o fluxo através do dente 
e meias cavas vizinhas e a força magnetomotriz correspondente. Esta relação pode ser 
determinada experimentalmente, por processo análogo ao adoptado no Cap. II; 

b) Considerar o campo associado às correntes do induzido na sua vizinhança imediata. Este 
campo pode, com razoável aproximação, ser tratado como uma perturbação ao campo ante- 
riormente determinado, considerando, para efeitos desta perturbação, o ferro como um meio 
anisótropo cuja permeabilidade na direcção do campo «não perturbado» e prêviamente 
conhecido seja y; , e cuja permeabilidade em direcções ortogonais seja p. 


5.2 Campo na zona do entreferro, dentes e cavas associado a correntes afastadas do domínio em estudo 
5.2.1 Método utilizado 


Para estudo do campo na zona do entreferro, dentes e cavas, seguiu-se processo análogo ao 
descrito no Cap. II, utilizando-se um modelo misto de papel condutor e rede de resistências. 

Fizeram-se ensaios correspondentes a diversos estados de saturação, por forma a poder inter- 
polar os resultados com aproximação satisfatória. 
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Verificou-se que a relutância do conjunto entreferro, dente e cava, é sensivelmente inferior ao 
valor obtido pelo processo corrente de cálculo, sendo a razão básica da discordância a diminuição 
acentuada de fluxo na região terminal do dente, a qual não é convenientemente tida em conta no 
processo habitual de cálculo. 

Considerou-se para f (H, B) a relação traduzida pelo gráfico da fig. 5. 


(1 
Fizeram-se ensaios correspondentes aos seguintes valores de Po ) e Bim J 
base 
ie base na gauss 
DS E neniá 

0,0000 (2) e 

0,0089 11 170 

0,0213 I2 600 

0,0393 13 185 

0,0565 13 830 


Para o estudo do caso limite de permeabilidade infinita utilizou-se um modelo em analogia 
inversa, analogamente ao exposto no Cap. II, em 5, tendo-se obtido o campo eléctrico representado 
por algumas equipotenciais na fig. 97. 


0.2.2 Resultados obtidos relativos à relutância do conjunto entreferro, dentes e cavas 


Os resultados obtidos foram referidos à relação 1! (3), definida no Cap. II, em 3 e 4,e são 
indicados no gráfico da fig. 98. Neste gráfico indicam-se também os valores obtidos pelo processo 
habitual de cálculo. Como seria de prever este conduz a valores por excesso. No gráfico da fig. 99 
indica-se o erro À (4) do processo usual de cálculo, referido ao valor obtido experimentalmente, 


o 


em função de (e) 


/ base 


Em valor relativo este erro é máximo para 


a 
246) Es - 0,02. 


base 


Este facto resulta de, nas condições vizinhas de máximo, a região do topo estar bastante mais 
saturada que a base. Dado que é na região do topo que o erro do processo usual de cálculo é 
bastante importante, o erro da relutância, no que respeita ao «campo no dente», tem valor relativo 
muito elevado. Para condições correspondentes a uma saturação mais acentuada, a variação de H ao 
longo do dente, em valor relativo, é menos pronunciada, sendo, nessas condições, menor o erro 
relativo global do processo usual de cálculo. 


(1) Usamos a notação definida no Cap. II. Assim (1) é a relação entre a permeabilidade magnética do 
wu / base 

vácuo e a permeabilidade média do dente na base, e Btm é o valor médio da componente radial da indução magnética 
ao nível do topo dos dentes (nos dentes e cavas). 

(2) Caso limite ideal de a permeabilidade do ferro ser infinita. 

(3) Relação entre a relutância do conjunto entreferro, dentes, cavas e a de um entreferro uniforme igual ao entre- 
ferro da máquina. 

(4) Definido no Cap. II, em 5 (erro do processo usual de cálculo em relação ao valor obtido experimentalmente, 
expresso em percentagem deste). 
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Além disso, de certa maneira, o facto apontado está relacionado com o facto de 


a relação f (H,B) considerada, apresentar um máximo para 


Ora, nas condições (em ordem de grandeza) em que se verifica o máximo de À, tem-se que, 


“A 


é da ordem de 0,03. 


E) 


do topo do dente, o valor de ( 


na região 
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5.2.3 Variação do fluxo magnético, do campo magnético e da indução magnética ao longo do dente 


de fluxo ao longo do 


dente, são representados na fig. 100, em que se indica a variação da relação entre o fluxo em diversas 


à variação 


+ 


Os resultados obtidos experimentalmente, no que respeita 


secções 


definidas pelo correspondente valor de 7 (1), e o fluxo na base 


) base 


do dente, ortogonais ao eixo, 


H 


Ho 


do dente, em função de ( 


Para 


algumas das curvas referidas apresentam um mínimo relativo. Tal facto resulta de, aproximadamente 


para essas condições, ser maior a variação relativa da permeabilidade ao longo do dente. 


A partir das curvas referidas pode-se calcular facilmente a variação do campo B ao longo do 


dente por forma a determinar as perdas magnéticas. Verificou-se que a variação de B em secções 


transversais é bastante pequena, sendo desprezável para efeitos da aproximação com que interessa 


calcular as perdas magnéticas. 


(1) Relação entre a distância da secção considerada à base do dente e a altura do dente. 
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5.3 Campo na zona do entreferro, dentes e cavas, na vizinhança das correntes a ele associadas 


5.3.1 Método utilizado 


Para estudo do campo na vizinhança imediata das correntes a ele associadas, considerou-se 
a analogia inversa, mais cómoda para a representação de um campo rotacional. 

Tipicamente este campo difere significativamente do anterior na cava em que se encontram os 
condutores respectivos e nos dois dentes vizinhos. 

“Considere-se, por exemplo, o caso representado na fig. 101 do campo associado a correntes 
na camada À y da cava a. Nas cavas b, c,.,., este campo tem configuração similar ao campo asso- 
ciado, por exemplo, a correntes no indutor. 

Seja o caso limite do ferro com permeabilidade infinita. O campo no entreferro no eixo dos 
dentes ab e bc, por exemplo, é muito aproximadamente radial. Consideremos o campo na cava e no 
entreferro entre os eixos dos dentes ab e bc. Dado que as condições fronteira definem univocamente 
o campo, o campo na cava b seria muito aproximadamente simétrico em relação ao eixo da cava. 

No caso da saturação a simetria é menos aproximada. Dado que, porém, o fluxo nos dentes 
ab e bc é da mesma ordem de grandeza, a distribuição do campo na cava b não se afasta muito de uma 
distribuição simétrica, coincidindo uma linha de força de H, aproximadamente, com o eixo da cava b, 

Na cava a e nos dentes vizinhos a distribuição do campo é bastante diversa da distribuição no 
resto da máquina. Na hipótese de uma distribuição de corrente nesta cava simétrica em relação ao 
seu eixo, o campo é antisimétrico em relação ao eixo. 

Supôs-se, para simplificar o estudo, que o campo no fundo da cava, na base do dente e na 
superficie do polo, era ortogonal a essas superfícies, hipótese que, estando o polo na vizinhança dos 
dentes prâticamente não saturados, é válida. Se tal não acontecer, porém, nada obsta a que as zonas 
do indutor e do núcleo do induzido vizinhas do campo em estudo sejam consideradas no modelo 
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Fig. 101 


experimental. Será apenas necessário um modelo um pouco mais extenso, e de que não nos foi pos- 
sível dispor. A justificação destas hipóteses é semelhante à apresentada no Cap. II, em 3. 

Supôs-se que o eixo das cavas vizinhas da cava em que se consideraram as correntes coincide 
com uma linha de força, hipótese que, pelos motivos acima expostos é razoavelmente aproximada. 
Aliás, a legitimidade de, com aproximação satisfatória, supor que os eixos das cavas vizinhas da- 
quela em que se consideram as correntes coincidem com linhas de força, foi verificada em face de 
resultados experimentais, conforme se verá mais adeante. 

Supôs-se que o campo no eixo da cava em que se encontravam as correntes era ortogonal ao 
eixo, o que, para o caso em estudo, é válido, desde que a distribuição das correntes na cava seja 
simétrica em relação ao seu eixo, o que, aliás, é usual. 


Assim sendo, tendo em conta as condições de simetria do campo, o seu estudo pode ser feito 
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a partir dum modelo que figure metade da cava percorrida pelas correntes em jogo, um dente e 
metade da cava vizinha, além do entreferro delimitado pelos dois eixos das cavas (fig. 102). 


cava 


— GINO Ss 
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e 
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Gsscociadoas 
ao campo 


Fig. 102 


Na hipótese de se não verificarem as condições de simetria indicadas o método exposto é tam- 
bém válido, com a necessidade, porém, de usar um modelo mais extenso. 

Usou-se uma representação baseada na analogia inversa. 

Dado que no entreferro e topos das cavas e dentes o campo é menos aproximadamente uni- 
forme representaram-se estas zonas por quadriculados correspondentes a menores dimensões con- 
forme, se indica ná fig. 103. 

Para o estudo deste campo utilizou-se uma rede de resistências (descrita no Cap. I, em 2.6, fot. 2). 

Tendo em conta os valores das resistências correntes no mercado e as suas características, e a 
montagem e processo de medida utilizados, decidiu-se considerar para resistividade plana corres- 
pondente à permeabilidade do ar 1 000 2. 

Nestas condições, as resistências são obtidas com relativa facilidade a partir de resistências 
correntes no mercado e são bastante estáveis. 

Fez-se o estudo para as hipóteses de o dente não estar saturado (permeabilidade práticamente, 
infinita) e dum estado de saturação correspondente a uma indução na base do dente da ordem de 
21 280 gauss (indução máxima no dente da ordem de 23 440 gauss). 

Em geral a distribuição de correntes na cava não é uniforme, quer em virtude de poder haver 
numa mesma cava correntes relativas a fases diferentes, quer devido ao efeito pelicular. Para o estudo 
deste efeito, convém conhecer, fundamentalmente, o campo associado a uma camada de correntes 
paralela ao fundo da cava, de pequena espessura, de distância ao fundo da cava arbitrária. Inte- 
ressa, portanto, estudar o campo associado a camadas de correntes de pequena espessura, nas con- 
dições atrás referidas. 

Os ensaios foram feitos admitindo que entre os condutores e o dente havia uma camada não 


, ' 1 
condutora (correspondente ao isolamento das bobinas) com E da largura da cava, de cada lado, 


sendo, numa secção da cava normal ao seu eixo, uniforme a distribuição de correntes. 
Utilizou-se a montagem indicada na fot. 7 e o esquema de ligações representado na fig. 104. 
As resistências da rede foram medidas e ajustadas com precisão, por forma que não diferissem 
mais de cerca de 0,5 "/y a 1 º/o do valor calculado. 
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As resistências R eram de precisão, bastante estáveis e de valores relativamente elevados, 
tendo sido medidas com rigor. 

Por intermédio de uma ficha os pontos A”, B', €”, D' e E” eram ligados aos nós da cava em 
que se pretendia injectar as correntes. 


CEAR TA TTI FEEL ATA 
ATI TAS II EIA NETS 
SEP PESERRARSSANEAREEEE 
= III TATI a 
==) JIIInDIITIII TITS 
== = SCE MBA RERASRA 
ESC raBELIBIienEssERERE 
ESSE cp SESBISGGSSASEEEE 
== =28 5 pa 
4 SAN III III II TIA 
=p pa BIESaanaRaEEa 
=SeEr mM RrairozanreERERE 
-EE ADA RIPAR RE E 
=EF Maas ERES naEnE 
=== CR mM RE Pa 
== SEE psEBISRESDRSREES 
= 
+++ 
Pr === = EP Dana BISaREna nana 
ESESAN JIIIILILIL III IT! 


Oscilos- 
8) -. cópio 


—s Resistência 


A resistência R,, de precisão, tinha um valor conhecido com rigor (1), por forma a permitir 
uma conversão numérica fácil da diferença de potencial entre dois pontos do modelo na relutância 
do correspondente tubo de força do campo magnético. 

As resistências R” eram ajustadas por forma que o potencial nos pontos A, B, C,D, E fosse 
o mesmo, impondo-se assim a condição de as correntes injectadas nos nós estarem na relação 
inversa das resistências R correspondentes. Por as resistências R terem um valor relativamente 
elevado este ajuste era bastante rápido. 

- A . . pa 1 — 

A resistência Ra era ajustada por forma que a tensão entre P; e P3 fosse 0 da tensão entre 

os terminais extremos do potenciómetro. Obtinha-se assim uma conversão fácil entre a diferença de 


CE e mm o me 


(1) Nos modelos correspondentes a permeabilidade do ferro infinita Ro, tinha um valor igual ao dobro da 
resistividade plana tomada para unidade (correspondente à permeabilidade do vácuo), e nos modelos correspon- 
dentes a uma forte saturação um valor igual a metade da resistividade plana tomada para unidade. 
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potencial entre dois pontos do modelo e a relutância do tubo de força correspondente, no campo 
magnético. 

Devido aos efeitos capacitivos parasitas, não foi possível utilizar corrente alternada, tendo-se 
empregado corrente contínua. Entre Pj e Ps aplicou-se uma tensão da ordem de 20 V. 


oo o 


' re e e e eg mg 
de as 


Fot. 7 


Convindo utilizar para elemento de zero um osciloscópio, recorreu-se a um interruptor síncrono. 
A tensão cujo anulamento se pretendia verificar era aplicada a um condensador. Por intermédio de 
um interruptor síncrono (50 c's) essa tensão era ligada periodicamente aos terminais de um oscilos- 
cópio (componente vertical do movimento do feixe). No osciloscópio obtinha-se uma função com 
a forma de impulsos rectangulares periódicos de amplitude igual à tensão cujo anulamento se pre- 
tendia verificar. 


5.3.2 Permeabilidade do dente práticamente infinita 


Para esta hipótese a zona representada na rede restringe-se ao entreferro e cavas. 

Fizeram-se ensaios correspondentes a correntes a diversos níveis da cava. Seja 7 a distância da 
camada ao fundo da cava, tomando como unidade a altura da cava. Os ensaios efectuados referiam-se 
a correntes nas camadas definidas pelos seguintes valores de 7: 0,05; 0,2; 0,4; 0,6; 0,8; 0,9875. 

Utilizou-se uma resistência R, == 2000 “ (fig. 104) de valor igual ao dobro da resistência 
plana unitária (1) (1000 4). 

Para cada uma das hipóteses consideradas mediu-se o potencial nos nós da rede, determi- 
nando-se, por interpolação, os pontos correspondentes a valores arbitrados de potencial. Unindo 


(1) Correspondente à permeabilidade do vácuo. 
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esses pontos obtiveram-se as equipotenciais representadas nas figs. 105, 106, 107, 108, 109 e 110, em 
que as equipotenciais são referenciadas pela correspondente leitura no potenciómetro (0 a 10º). 

Qualitativamente os resultados obtidos são facilmente justificáveis. As linhas de força do campo 
magnético (equipotenciais do modelo) afastam-se tanto mais da ortogonalidade ao eixo da cava 
quanto mais próximas do entreferro. A ortogonalidade é razoavelmente aproximada excepto na vizi- 
nhança do entreferro. 

Para verificação da legitimidade de ter suposto que os eixos das cavas vizinhas daquelas em 
que se encontravam as correntes coincidia com linhas de força do campo magnético, representou-se 
na rede de resistências a hipótose limite de simetria, impondo que o eixo de um dente e o eixo da 
cava coincidissem com linhas de força, tendo-se obtido o campo representado na fig. 111. Este campo 
coincide muito aproximadamente, na zona da cava que representa, com os campos rspresentados nas 
nas figs. 105 a 110 (e com o representado na fig. 97), o que confirma a legitimidade da hipótese feita. 

Sejam dois pontos P; e P; do modelo, de potenciais V' (Q1) e V' (Q:) tomando como unitária 
a tensão entre os terminais extremos do potenciómetro (nas condições experimentais atrás descritas), 
e sejam: 


P,;z o fluxo em maxwell do tubo de força do vector indução magnética definido por Pj e Ps 
e por dois planos ortogonais ao eixo da máquina e à distância |, em centímetro ; 


fm o fluxo do vector densidade de corrente, em ampere, na camada de corrente associada 
ao campo magnético, na metade da cava representada no modelo ; 


| a extensão do campo perpendicularmente à secção representada no modelo, em centímetro. 
Seja ainda 


249) V(Q)=V'(Q) — 0,1000. 


Tem-se 


250) MM E = 82 [(V(QU—V (QI =8 7 [V(QU --V(QUI. 


m e 


Fundamentalmente, para uma análise quantitativa, interessa conhecer o fluxo médio asso- 
ciado a uma camada definida por um valor 74 de 7, fluxo esse associado a correntes numa camada 
definida por um valor 72 de 7, para quaisquer valores de 7 compreendidos entre O e 1. 

Esse problema reduz-se à determinação do valor médio Vm (71, 73) do potencial numa camada 
definida por 7==t4, sendo esse potencial associado a correntes na camada definida por +=". 

Para a determinação deste valor médio utilizaram-se os valores obtidos nos ensaios para os 
potenciais nos nós da rede, com factor de ponderação 50 º/o no nó imediatamente vizinho do 
dente, para ter em conta a espessura de isolamento entre os condutores e o dente. 

Os valores de Vm (t1, 72) assim calculados foram marcados num diagrama, de coordenadas 
Tt, 72, determinando-se por interpolação as curvas correspondentes a um mesmo valor de 
Vm (r1, 7a) (fig. 112). 

A partir, por exemplo, deste diagrama, e tendo em conta as restantes condições da máquina, 
é fácil calcular a distribuição de correntes na cava, e, conhecida esta, a distribuição resultante do 
fluxo, e portanto as reactâncias associadas ao fluxo nas cavas. 

Deve notar-se que o método experimental utilizado é igualmente aplicável à hipótese de se 
não admitir a uniformidade da distribuição de correntes em secções ortogonais ao eixo, permitindo 
até determinar essa distribuição. Há, apenas, que considerar, em cada um dos ensaios experimentais, 
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Fig. 106 


o — Nós em que as correntes foram introduzidas no modelo. 
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Fig. 107 
() — Nós em que as correntes foram introduzidas no modelo. 
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Fig. 108 


correntes foram introduzidas no modelo. 
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* PerbunanN - resistente ao calor e ao envelhecimento 


Artigos de borracha feitos da borracha butadieno-acrilnitrilica Perbunan N, da Bayer, são muito resistentes ao calor e 
ao envelhecimento, resistem muito bem ao ataque de óleos, gorduras e gasolina, e são também resistentes ao desgaste. 


Os nossos tipos de Perbunan N são excelentemente apropriados para a fabricação de artigos técnicos moldados, 
resistentes aos óleos e à gasolina, especialmente para solas e tacões resistentes ao desgaste, assim como para 
cintas de travões, para rolos, borracha de carimbos e telas de cauchu, resistentes aos óleos, para correias 
trapezoidais, correias transportadoras, tecidos gomados, bainhas para cabos eléctricos e mangueiras, especialmente 
mangueiras de gasolina. 


Produtos de alta qualidade exigem um controle permanente. O Perbunan N é controlado continuamente para 
garantir uma alta qualidade sempre igual do material. 


PERBUNAN N a Borracha Nitrílica da BAYER 


Representante em Portugal: 

S.A.R.L. QUIMICOR 

Lisboa: Rua Sociedade Farmacôutica, 3 
Tel. 555496 e 555893 

Porto: Rua de Santos Pousada, 441 
Tel. 54141 
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MAKROLON - 

o Policarbonato 
da BAYER 
Qualidade e beleza 
de formas 


Três propriedades, a de ser inquebrável, 

a insipidez absoluta e a de ser indeformável em 
água fervente, são a razão porque, hoje 

em dia, se fabricam de Makrolon, 

o policarbonato da Bayer, objectos de utilidade 
doméstica, de toda a espécie e em grande 
número. As peças moldadas deste termo- 
plástico são resistentes a altas temperaturas - 
até 135º C. - e repelem os corantes. 

As substâncias corantes contidas no café, 

no chá, nos sumos das frutas e até mesmo nos 
«batons » não aderem ao Makrolon. 


Em todo o mundo se empregam os «biberons » 
de:Makrolon. Louças de qualidade, moldadas 
de Makrolon aliam o belo ao útil. 


Na Feira dos Plásticos em Duesseldorf, de 12 a 
20 de Outubro de 1963, damos uma visão da 
nossa capacidade de produção. Pode falar 


connosco sobre plásticos, no Pavilhão E / 
rés-do-chão, assim como no stande ao ar livre, S.A.R.L. QUIMICOR, Representante em Portugal 


em frente do Pavilhão E. No Pavilhão F 1 Lisboa: Rua Sociedade Farmacêutica, 3, Tel. 42194 
damos informações sobre borracha sintética. Porto: Rua de Santos Pousada, 441, Tel. 54141 
Teremos muito prazer em receber a sua visita, 

que gostosamente aguardamos. 


| : 
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no modelo. 


Nós em que as correntes foram introduzidas 


O 
TRONIOA 


728 


oo: 00014 


0:01 
pros Foo: 


osso: 


+ 
oo 
ooz! | 
00Fr 


c08: * 


TRONIOS 
729 


Fig. 111 
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em vez de uma distribuição de correntes numa camada ortogonal ao eixo, uma corrente num ponto 
definido por 7 e pela sua distância é ao eixo. O ponto (filete ortogonal ao plano considerado) a que 
se refere o fluxo (ou V), e o ponto correspondente às correntes associadas ao campo, são definidos 
por 7 e pela sua distância ao eixo da cava. 

Os resultados dos ensaios experimentais permitem uma melhor aproximação que a correspon- 
dente a considerar apenas Vm (71, 73), considerando fluxos associados a distribuições uniformes de 
correntes, mas definindo os pontos a que se refere o fluxo pelas coordenadas 7 e £. Pode-se, assim, 
estabelecer uma «melhor» aproximação para a distribuição de correntes, correspondente à calculada 
«a partir» de uma distribuição uniforme de correntes em secções ortogonais ao eixo. Como é óbvio 
trata-se de uma correcção de segunda ordem em relação à não uniformidade de corrente ao longo 
do eixo. 

5.3.3 — Dente fortemente saturado 


Como atrás se disse, o campo na vizinhança das correntes a ele associadas, na hipótese de 
saturação do ferro, foi considerado como uma perturbação ao campo associado às restantes corren- 
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ão entre a permeabilidade e a permeabilidade diferencial tradu- 


turbado e bem assim a direcção do campo. 
113. 
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tes na máquina, admitindo que, para efeitos desta perturbação, o ferro se comporta como um meio 


anisótropo de permeabilidade, na direcção do campo não perturbado, pj, e em direcções ortogonais, £. 


Houve assim que obter prêéviamente as condições correspondentes ao campo não perturbado, 


para o que se estudou, pelo método descrito no Cap. Il em 3.1, uma distribuição do campo, corres- 


pondente a uma indução na base de dente da ordem de 21 280 gauss. 


A partir da distribuição assim obtida, ficou-se a conhecer a distribuição de permeabilidades 


no campo não per 


Verificou-se não haver êrro sensível em considerar para direcções principais do tensor per- 


meabilidade correspondente à perturbação as direcções radial e tangencial. 


A partir da distribuição de permeabilidades obtida, calculou-se a distribuição de permeabili- 


dades diferenciais a partir da relaç 


zida pela curva da fig. 


Calcularam-se as resistências do modelo, em analogia inversa (pelo método indicado no Cap. I, 
em 1 e 2), correspondendo as resistências radiais à permeabilidade  (resistividade proporcional ao 


inverso da permeabilidade) e as resistências tangenciais à permeabilidade diferencial pi. 


No gráfico da fig. 114 indica-se a variação da permeabilidade » e da permeabilidade diferencial 


Hi, em média, ao longo do eixo. 
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Fig. 114 


Procedeu-se então ao estudo do campo associado a correntes na cava em camadas ortogonais 


pelo método atrás descrito, tomando para R, (fig. 104) o valor 500 2), igual a metade da 


resistividade tomada para unitária (1), tendo-se feito ensaios para correntes correspondentes aos 
mesmos valores de 7 que na hipótese de permeabilidade do dente infinita. Os resultados obtidos são 


representados nas figs. 115, 116, 117, 118, 119 e 120. 


ao eixo, 


Qualitativamente, estes resultados têm uma interpretação fácil. É notória, no entanto, a dife- 


rença entre este campo e o campo na hipótese de não haver saturação. 


Para verificação da legitimidade de admitir que o eixo das cavas vizinhas aquela em que se 
encontram as correntes coincide com uma linha de força do campo magnético, fez-se o estudo do 


à permeabilidade do vácuo. 


(1) Correspondente 
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Fig. 115 


() — Nós em que as correntes toram introduzidas no modelo. 
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Fig. 116 


() — Nós em que as correntes foram introduzidas no modelo. 
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Fig. 117 


(O — Nós em que as correntes foram introduzidas no modelo 
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Fig. 118 


(O — Nós em que as correntes foram introduzidas no modelo, 
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Fig. 119 


(O — Nós em que as correntes foram introduzidas no modelo. 
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Fig. 120 


(o — Nós em que as correntes foram introduzidas no modelo. 
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Fig. 121 


campo «de perturbação» na hipótese de simetria, nas mesmas condições de saturação, impondo que 
o eixo de um dente e o eixo da cava vizinha coincidissem com linhas de força, tendo-se obtido o 


campo representado na fig. 121. Da comparação desta figura com as figs. 115 a 120 ressalta a 


razoável legitimidade da hipótese feita. 


MIA NAL ALINALNAL NI TS. 


ea NAN NV IN INN NAL SS 


A NAINNAM NNININCEO 


ISSN NS AN NINE NO 


RL NSNEANE NENE 
INNER NT 


TISKEL KEN. NENE 
SNIS: 
TONENCENCN NE DEITDSEE 
SNS SNSNESHEO 


ISNTNT NT NT NT RINS 
ENE PRENDE NREcE 
N Ele 
ANINALMIAL LA ALNATANT 


ARIAL RLIBIRARREESS o 
o a q Es o o » m eu q E ias 
= S RS q q S s S S Q * gs 


Tem-se, com notação idêntica à definida em 5.3.2, sendo V” (Q) o potencial de um ponto do 


modelo (tomando como unitária a tensão entre os terminais extremos do potenciómetro), 


=s8r| V (O — (09 |, 


27 | V(Q)—V” (Og | 


sendo 


| V“ (Q) — 0,100 | 


a 
4 


252) V(Q) 
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Também por processo análogo ao descrito em 5.3.2 obteve-se a função Vm (71, 73) representada 
no diagrama da fig. 122. 
É interessante a simetria do diagrama obtido em relação à recta de equação 


253) tÚ=u. 


Essa simetria corresponde à igualdade dos coeficientes de indução mútua entre dois circuitos, 
em regime linear, e está implícita na representação por intermédio de uma rede de resistências 
merecendo ser analisada com algum detalhe. 

Considere-se uma rede de resistências (elementos lineares) em corrente contínua (em regime 
estacionário) com m nós distintos (ligados por resistências não nulas) e referenciados por números 
de ordem 0,1,2,...m—1. 

Considere-se um desses nós, referenciado pelo número de ordem O, em relação ao qual sejam 
referidos os potenciais dos restantes. 

Seja i, a corrente introduzida no nó genérico «, de potencial V,. Seja Rg, a resistência que 


liga dois nós genéricos P e 7. 
A relação entre as correntes i, e as tensões V, pode ser expressa pelo seguinte sistema 
de equações lineares : 


i=CuV + CaVa +... + Cim- Yna 
is = Co Vi + CuVa +... + Com Vma 


da == Cia Vi + Cu, Va + ... + C, (m-1) Vm- 


254) 
im =C(m-y Yi + Ca Va +... + Com my Vm-t 
E e eis De q ESB) 
Vo = 0 
ou 
m—l 
= >, Cu V3 
B=1 
254º + ca 
lo =>, iz 
Vo= O 


Para calcular os coeficientes C,; consideremos um valor particular de É, arbitrário e admita- 
mos a hipótese de ser 
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Tem-se 


1 
257) la — p= Ci = 5 R 


258) ictp=Cp=— 
a+B Rag 


Tem-se, portanto, que 


259)  Cg==Cha - 


Consideremos a matriz dos elementos C,g e o determinante correspondente À (simétrico). 
Seja 4, o produto por (— 1)*+É do menor de primeira classe desse determinante relativo ao 


elemento C,g. 
Tem-se 


m—l 
> 8, B la 
a=1 


a! 


260)  Va= 


relação que se pode pôr sob a forma 


m=1 
0) 4-5 p,,1 
[rÀ A 
gi==1 
sendo 
A 
0) Dig E = Doe 


Sejam dois conjuntos de nós (não compreendendo o nó 0), um de elementos pr ,..,Pq,- «Pr; 


e outro de elementos s,,..., St,..., Su. Sejam kpy,..., kpq 1 ++ -o a guess x a q ss j 
ks; ,...,y Ksy, parâmetros constantes. 
Tem-se 


262) Sky ku Do Dk, toa Dep 
q t t q 


Esta igualdade traduz, em analogia inversa, a condição de igualdade dos coeficientes de indu- 
ção mútua (1) entre dois circuitos, um correspondente a condutores ropresentados nos nós p e outro 
a condutores representados nos nós s, sendo o número de condutores de qualquer dos circuitos figu- 
rado nos diversos nós em geral diferente (proporcional aos coeficientes k). Os fluxos a que se refe- 
rem estes coeficientes de indução consideram-se definidos em relacão ao mesmo nó O. 


5.3.4 — Comparação dos resultados obtidos nas hipóteses de não saturação e de saturação do ferro 


Para comparar os resultados obtidos nas hipóteses de não saturação e de saturação do ferro, 
calculou-se, a partir dos valores experimentais, a função 


) 


263) = f (1), 


mt 


(1) Um fluxo do vector indução magnética através dum dado percurso é representado por uma circulação do 
campo eléctrico através do mesmo percurso no campo eléctrico análogo inverso. 
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sendo: 


P o valor médio do fluxo do vector indução magnética entre o eixo da cava vizinha daquela 
em que se consideram as correntes e os pontos de uma secção da cava ortogonal ao eixo 
(secção essa definida pelo correspondente valor de 7), na hipótese de uma distribuição uni- 
forme de correntes na zona da cava indicada a tracejado na fig. 123; 


Fig. 123 


fmt a corrente total na cava associada ao campo magnético. 

Obtiveram-se as funções representadas graficamente na fig. 124. 

Conclui-se, em primeiro lugar, destas curvas, a importantíssima influência da saturação no que 
respeita ao fluxo na vizinhança imediata das correntes (considerado como uma perturbação ao campo 
magnético associado a correntes afastadas do domínio em estudo). 

Para as condições a que se referiam os ensaios efectuados tem-se que, por exemplo, na hipó- 
tese considerada de distribuição de correntes: 


a) Os fluxos, definidos em relação aos eixos das cavas vizinhas, e relativos à zona da cava em 
que se encontram as correntes, nas duas hipóteses de saturação a que se referiam os ensaios, 
estão numa relação da ordem de 11:1 


b) Os fluxos entre o fundo da cava e o entreferro estão numa relação da ordem de 15:1. 


Estes resultados mostram que, na hipótese de saturação apreciável, não é admissível supor, 
para efeito dos fluxos na vizinhança das correntes, que o ferro se comporta como um meio de per- 
meabilidade infinita. 

A variação de fluxo ao longo da cava determina a distribuição de correntes ao longo da cava, 
a qual, na hipótese de se não usarem artifícios construtivos, é bastante diversa de uma distribuição 
uniforme. Nestas condições, na hipótese de saturação, o cálculo efectuado supondo a validade de 
considerar infinita a permeabilidade do ferro, conduz a resultados bastante pouco aproximados. 

Os diversos problemas relacionados com a distribuição de fluxos na cava (cálculo da distri- 
buição de correntes e das reactâncias) podem ser fácilmente resolvidos a partir do conhecimento de 
V (r4, 73), se necessário por aproximações sucessivas. 

Experimentalmente determinou-se V (t1 , 73) apenas para dois estados de saturação. 

Tendo em conta o andamento do campo (figs. 115 a 120), constitui, no entanto, aproximação 
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razoável (1), interpolar entre os valores obtidos nos dois ensaios, considerando para variável de inter- 
polação o valor médio do inverso da permeabilidade diferencial ao longo do dente. 
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Designando por Va (r1, 72), Vb (T1, 73), Vw (71, 73) os valores de V (m , 73) respectivamente 
nas condições ensaiadas, sem e com saturação, e noutro qualquer estado, tem-se 


264) Velu,T)=(1—Ê)Va(m,7.) +EVo(m,T), 


sendo É um parâmetro dependente do estado de saturação e proporcional ao valor médio do inverso 
da permeabilidade diferencial ao longo do dente. 

Por forma a permitir utilizar com mais facilidade os resultados obtidos, indica-se o valor de É 
em função de Bim, valor médio da componente radial do vector indução magnética do campo «não 
perturbado» (sem ter em conta as correntes na cava em estudo), ao nível do topo dos dentes 


(fig. 125). 


—— —— e 


(1) Pelo menos à falta de mais resultados experimentais. 
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Na fig. 126 representa-se, para diversos valores de Bim, a função 


p 
265) ————— =t (1), 
) a (7) 


na hipótese de uma distribuição uniforme de correntes na parte da cava indicada a tracejado na 
fig. 123. 
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6 Possibilidade de obtenção rápida de resultados experimentais 


Os valores experimentais obtidos cobrem uma gama relativamente reduzida dos diversos pará- 
metros, podendo surgir deles a ideia da morosidade do método. Tal resultou, no entanto, da carên- 
cia de meios experimentais. Na medida em que se disponha de equipamento experimental adequado, 
os ensaios podem tornar-se bastante mais rápidos, além de que, em estudos de carácter sistemático, 
é possível uma muito sensível diminuição do tempo necessário à obtenção dos valores mais impor- 
tantes, com uma cuidada programação. 

Vejamos um exemplo: 


Na medida em que, no método experimental descrito em 5.3, interesse unicamente conhecer 
Vm (71,72) com factores de ponderação conhecidos para os diversos nós correspondentes a 74, basta 
determinar o valor médio do potencial nos nós correspondentes a 71, não interessando os potenciais 
de cada nó. 


Esse valor médio pode ser obtido com uma única leitura utilizando o esquema da fig. 127, em 
que 1,2,...,i,... são os nós da rede cujo potencial médio, com factores de ponderação conheci- 
dos, se pretenda determinar, e R4,..., Rj,..., são resistências auxiliares inversamente proporcio- 
nais aos factores de ponderação relativos aos nós a que são ligadas. 
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Sendo nula a corrente i,, tem-se 


Se as resistências R; forem muito superiores às resistências entre os nós da rede, os potenciais 


la leitura 


ões, pe 


ivelmente alterados pela ligação das resistências Rj. Nestas condiç 


do valor de V, tem-se directamente o valor médio pretendido. 


ão são sens 


— 


V, n 


7 Obtenção dos valores globais a partir do conhecimento do comportamento local 


das máquinas 


7.1 Introdução 


Consideremos por exemplo um alternador trifásico. À parte alguns factores a que nos referi- 
remos mais adiante, nada obstando aliás, salvo a maior complexidade analítica, a que sejam também 
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tidos em conta, o comportamento electromagnético desse alternador pode ser obtido a partir das 
seguintes relações : 


O, = (6, lo sM,di 38) 


267) O =f(a,io,iMu,is,is) 
da=f(z,io,iu,iz, is) 


| Da = E (4, adt, O, DB), 


| tnalrumento de zera 
sendo Fig. 127 

b, o fluxo através do enrolamento indutor ; 

D,, Pa, P; os fluxos através dos enrolamentos das fases 1, 2 e 3; 

a um parâmetro definidor da posição angular do rotor; 

io a corrente no rolamento indutor; 

Hu, iz, is, as correntes nos enrolamentos 1, 2 e 3 do induzido. 


Nas expressões atrás indicadas supôs-se, para simplicidade de exposição, que os quatro enrola- 
mentos da máquina podiam ser tratados como circuitos bem definidos (supondo, por exemplo, uma 
distribuição uniforme de corrente na secção dos condutores), que não havia enrolamento amortecedor, 
que se podia desprezar, para efeitos de valores totais de fluxos, as correntes no núcleo magnético e 
a histerese do ferro. Trata-se de simplificações legítimas para uma exposição global; na devida altura 
deverão ser analisados os erros correspondentes, sendo ou não abandonadas essas hipóteses. 

O comportamento dinâmico da máquina e as características dos circuitos a que esteja ligada, 
fornecem um conjunto de relações entre as variáveis consideradas, conjunto que, simultâneamente 
com as relações (267), permite definir o comportamento da máquina. 

Não se trata, evidentemente, de um problema em geral simples. Trata-se, porém, de um pro- 
blema diferente do que abordámos (embora com ele relacionado), que é o comportamento local das 
máquinas eléctricas. 

Para exemplificar como poderão ser tidos em conta os resultados obtidos pelo método que expu- 
semos, vejamos como podem ser determinadas as relações (267) referentes a uma parte da máquina 
considerada plana, e em que seja legítimo desprezar a saturação do ferro nos polos indutores e no 
núcleo do induzido. 


7.2 Fluxo através do enrolamento indutor 


O fluxo ?, através do enrolamento indutor, na zona da máquina tratada como plana (de espes- 
sura |, em centímetro), pode, em geral, ser obtido pela soma de duas parcelas: 


a) O fluxo 9, associado à zona debaixo dos polos. Por exemplo, na hipótese de um alternador 
de polos salientes, esse fluxo pode, em geral, ser determinado a partir da soma dos fluxos 
9';(1) nos dentes, definidos por exemplo entre os eixos das cavas vizinhas no fundo das cavas. 


Será 


268) P,=mn b> Ub', 
j 


sendo 


No o número de espiras do enrolamento indutor; 
U =U (x,j) uma função que tem o valor 1 para os valores de « e j tais que os dentes res- 
pectivos estejam debaixo do polo e o valor O para os valores de x e j tais 


(1) Representamos por j um índice que permita identificar os dentes da máquina. Para simplicidade de notação 
definimos j pela distância angular do eixo dos dentes a um ponto fixo do estator, tomando como unidade o passo dos 
dentes e convencionando um sentido positivo. 
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que os dentes respectivos estejam entre os polos; esta função U é facilmente 
determinável. 


O problema do cálculo de ?', reduz-se, portanto, ao cálculo de ';. 
b) O fluxo P”, associado à zona entre os polos, que é aditivo, e que pode ser obtido por uma 
expressão da forma 


H 


269) -  =No | fm E = > mo É” (De) | , 
s 


sendo 
fmm o fluxo do vector densidade de corrente associado às bobinas indutoras; 
Pp” 
É uma grandeza que traduza o valor de ——º- na hipótese de só haver correntes no 
mm - 
enrolamento indutor (por exemplo a grandeza analisada em 4.2.3); 
Eme o fluxo do vector densidade de corrente associado a uma cava genérica c; 


tt 


f' (Dc) uma função que traduza o valor de na hipótese de só haver corrente na 
me + 
cava definida por c, dependendo Dc!) de c e de 2. 


Tem-se, portanto, que = se pode reduzir à forma 


H 
p é 


270) = Aoilo + Arin+ Asis+ Asis, 
sendo As, Ar, A» e As funções de «2, 
O fluxo através do enrolamento indutor será 


Do ea dr d Po 


271 ” 
| | | 


7.3 Fluxo através dos enrolamentos do induzido 


Seja G(2) a distância angular de eixo de um dente ou de uma cava a um ponto fixo do esta- 
tor, tomando como unidade o passo dos dentes, e convencionando um sentido positivo. 

Seja o fluxo Pk (k==1, 2 ou 3) através de um enrolamento k, constituido por my espiras nas 
cavas definidas por valores ri e si de G, na nas cavas definidas por r2 e s2,..., ni nas cavas defini- 
nidas por r; e s;,..., sendo ri; si. 

Sejam : 


04"; o fluxo através do dente definido por j, entre os eixos das cavas vizinhas, no 
fundo das cavas; 


Qr-+o5 o fluxo médio entre os condutores da cava ri; pertencentes ao enrolamento 
em causa e o eixo da cava (no fundo da cava) definida pelo valor r,!-1 de G; 
i 


QD's —0,5 o fluxo médio entre o eixo da cava (no fundo da cava) definida pelo valor 
Ss; —1 de G e os condutores da cava si pertencentes ao enrolamento em causa, 


Tem-se 


si — 1,5 + 
272) Dj — 3 je ni [, 4 0,5 + E — 05 + > 0, | . 
Í 


i=r]+1,5 


— [000000 00 


(1) Sendo De , por exemplo, a distância angular do eixo da cava a um eixo entre polos, tomando como unidade o 


passo dos dentes, e convencionando um sentido positivo. 
(2) Os valores de G e D correspondentes a um mesmo ponto estão relacionados por intermédio de a, parâmetro 


definidor da posição relativa do rotor e do estator. 
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Portanto, o problema da determinação de Pk está reduzido ao da determinação de ?,, 'P',, +o,5 
e 'b's, — 0,5, em função de «, io, ir, is, is. 


7.4 Determinação dos fluxos b; 


' - mo dd ' 
Comecemos por considerar o caso de não haver saturação do ferro, sendo ?; o valor de 9; 
nesta hipótese. Tem-se 


tt 


bp 
273) =fmm É“ (Dj) + 5) fme (Dj, D;— Do), 


Cc 


sendo 
P; o fluxo, em maxwell, através do dente definido pelo valor ;j 
de G, entre o eixo das cavas vizinhas, no fundo da cava, 
e no domínio definido por | (em centímetro), na hipótese de 
não haver saturação, e associado às correntes no indutor; 
ff” (D;) ou f” (j, a) uma função que traduza o valor de Es na hipótese de só 


haver correntes nos enrolamentos indutores (por exemplo a 
função definida em 4.2.3) ; 
! 
fº (D;, DD) ou £* (j, 2, j-c) uma função que traduza o valor de 


na hipótese de só 
mc + 


haver correntes na cava definida pelo valor c de G (por 
exemplo a função definida em 4.4.5). 


t 


Tem-se, portanto, que = se pode pôr sob a forma 


t 


Pp 
274) = Di GDE G,0).izt E (j,0).is. 


e 
Sendo 2”, definido nas mesmas condições de P;, com a diferença, porém de se considerar a 
saturação, tem-se 


275) ae S; (A, ) , 


tt 


sendo 5; A a) a relação entre o fluxo nas hipóteses de saturação e não saturação. Esta relação, 


na zona debaixo dos polos, pode ser obtida com a maior facilidade a partir, por exemplo, da função 


276) T=T (Bm) 


definida no Cap. II, em 3.4. Na zona entre polos não afectada pela saturação, tem-se 


277) S5,=1. 
7.5 Determinação dos fluxos D+ 0,5 € Me 0,5 


Seja, com notação análoga à definida atrás (1) 


+ 


278) ni à + 0,5 | = Emm ii “ r pá + > Eme É'º (D; ' D; — De) ' 
4 — 0,5) | — 0,5 cr; 
ea | 


—e a 


(1) Consoante os casos se deverá considerar, respectivamente, o primeiro ou o segundo dos valores incluídos 
no (7. 
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sendo 'D! I ri 0,5 | o fluxo, em maxwell, através do dente definido pelo valor | H + a de G, entre 
sji — 0,5 *% — 0,5 


os eixos das cavas vizinhas no fundo das cavas, no domínio definido por |, fluxo esse associado a todas 


E o A “i cê é 
as correntes na máquina que não as na cava definida pelo valor | h de G (1), na hipótese de não 
à 


se considerar a saturação. 


Seja pr r +0,5|] o fluxo definido em condições análogas a 'P”, com a diferença de 


S; — 0,5 
considerar a saturação. 


Tem-se 


Sj — 0,5 


e” | "+ sa 
599) rj meo Ee | rj e É r, +0,5 | , a 


S; — 0,5 Si — 0,5 Si — 0,5 | 


Para obter 'b' [ri 40,5] há que sobrepor a “MP” [r; 0,5] O campo associado a correntes 
Sj — 0,5 Sj — 0,5) 


r 
na cava | k ) , por exemplo pelo método exposto em 4.4.6 e em 5.3. 


7.6 Cálculo do comportamento global das máqninas 


Uma vez conhecidos 9, Pr+os e P's, —o,5 em função de &, io, iu, iz e is, tem-se, substi- 
tuindo em (268), (271) e (272), um conjunto de expressões da forma (267). Deve notar-se que, sendo 
o alternador de construção simétrica, as relações (267) (1) satisfazem a diversas condições de simetria 
relativas às variáveis «, i, iz e iz, sobremaneira evidentes para nos dispensarmos da sua justificação 
e análise. 

Como se disse atrás, as relações (267), conjuntamente com as que traduzem o comportamento 
dinâmico da máquina e as relações entre as diversas grandezas e o tempo impostas pelos circuitos 
a que esteja ligada, permitem definir o comportamento global da máquina. 

Assim, por exemplo em regime permanente, igualando « ao desvio angular em radiano entre 
um ponto do rotor e um ponto do estator, e sendo 


p o número de pares de polos do alternador, 

f a frequência (eléctrica) do rotor, 

0= 27% f 

Vo a tensão entre os terminais do circuito indutor, 
To a resistência do circuito indutor, 


Vi, Vs e v; as tensões entre os terminais dos enrolamentos do induzido, 
ry, Te Tr, as resistências dos enrolamentos do induzido, 


ter-se-á: 


280) aÀ =— nt+a, 
p 
281) 4, = const. 


(1) Supõe-se que os condutores da fase em relação à qual se calcula o fluxo têm uma distribuição simétrica em 
relação ao eixo da cava. O método exposto é fácilmente generalizável à hipótese de tal se não verificar. 
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eportagem General Electric sobre energia n 
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O primeiro reactor do mundo de vapor sobreaquecido, poderia 
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reduzir o custo da energia nuclear 


(Em cima) O reactor experimento! de vapor sobreaquecido ESADA em 
Vallecitos (em primeiro plano), Pleasanton, Califórnia, E. U. A. que foi 
construído pela General Electric e pela Empire State Atomic Development 
Associates, Inc, precursor de grondes instaloções de vapor sobreaque- 


cido em centrais nucleares, que estorão prontas por volta de 1968. 


A direita) A central nuclear de Garigliano (150.000 kw ), situada o 

milhas ao Norte de Nápoles, na ltólio, pertencente à Societô 
Elettronucleare Nazionale" (SENN), depois de instalada q linha de 
transporte de energio eléctrico; esto central estó já abastecido de 
combustível nuclear. O reactor de água fervente de ciclo duplo 
começou a funcionar em Maio. 


Estes reactores General Electric-— um deles uma instalação 
experimental para experiências 'de combustiveis onde cien- 
tistas americanos investigarão o uso do vapor sobreaquecido 
em centrais nucleares como processo de reduzir o custo da 
energia; o oulro a maior instalação produtora de energia, 
no seu género, na Europa — são apenas dois dos cinco reac- 
tores GE que entrarão em funcionamento este ano. 

Os outros três são: o Reactor Experimental Japonês de 12,5MW, 
em Tokai-Mura (carga de combustível prevista para Julho); 
o de 50/60MW, da Central Nuclear de Humboldt Bay, na 
Colifórnia (atingiu o ponto critico em Fevereiro) e o de 
50/75MW, da Central Nuclear de Big Rock (atingiu a potência 
máxima em Março). 

À princípio, o reactor experimental de vapor sobreaquecido 
ESADA funcionará com a potência 12.500 kW. Será utilizado 
pelos investigadores para experiências com combustíveis e 
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